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We investigate the transient dynamics of a lumped-element oscillator based on a dc supercon¬ 
ducting quantum interference device (SQUID). The SQUID is shunted with a capacitor forming 
a nonlinear oscillator with resonance frequency in the range of several GHz. The resonance fre¬ 
quency is varied by tuning the Josephson inductance of the SQUID with on-chip flux lines. We 
report measurements of decaying oscillations in the time domain following a brief excitation with 
a microwave pulse. The nonlinearity of the SQUID oscillator is probed by observing the ringdown 
response for different excitation amplitudes while the SQUID potential is varied by adjusting the 
flux bias. Simulations are performed on a model circuit by numerically solving the corresponding 
Langevin equations incorporating the SQUID potential at the experimental temperature and using 
parameters obtained from separate measurements characterizing the SQUID oscillator. Simulations 
are in good agreement with the experimental observations of the ringdowns as a function of applied 
magnetic flux and pulse amplitude. We observe a crossover between the occurrence of ringdowns 
close to resonance and adiabatic following at larger detuning from the resonance. We also discuss 
the occurrence of phase jumps at large amplitude drive. Finally, we briefly outline prospects for 
a readout scheme for superconducting flux qubits based on the discrimination between ringdown 
signals for different levels of magnetic flux coupled to the SQUID. 


I. INTRODUCTION 

Superconducting circuits composed of Josephson junc¬ 
tions have been the subject of intense research for the 
past few decades for their importance in understanding 
the fundamental aspects of quantum mechanics as well 
as for their potential application towards quantum infor¬ 
mation processing and computing These /rm-sized 

devices have been shown to exhibit macroscopic quantum 
tunneling Si, quantized energy levels S and superpo¬ 
sition of states in a quantum bit (qubit) S-S]. A central 
^plication of Josephson devices, classical or quantum, 
[9| is measurement. There have been many advances in 
utilizing Josephson devices such as a dc superconducting 
quantum interference device (SQUID) for qubit readout 
[ifl-US- In these experiments, the SQUID forms part of 
a resonant oscillator circuit that is coupled to a qubit. 
Some of these Josephson devices configured as amplifiers 
hav e ap proached the quantum limit in noise performance 
[lJ-[^. All of these applications of Josephson junctions 
depend at some level on the nonlinearity of the junction 
response. 

In this paper we investigate the temporal dynamics of 
a nonlinear SQUID resonant circuit. Studies of nonlin¬ 
ear oscillator dynamics under continuous excitation have 
been done previously fl^ . [23 - [^ . and, in fact, many types 
of superconducting qubits, such as the transmon |28j or 
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the phase qubit @ , are nonlinear oscillators that are typ¬ 
ically driven with resonant pulses. However, we are not 
aware of any experimental or theoretical work to date 
on the transient dynamics of nonlinear oscillators under 
pulsed excitation in the time domain. Here we present 
time domain measurements of the decaying voltage os¬ 
cillations from the SQUID oscillator after a brief exci¬ 
tation. The SQUID potential, and hence the resonance 
frequency, can be tuned by changing either the bias flux 
or bias current; here we focus on the variation with re¬ 
spect to flux, while no dc bias current is applied. 


The paper is organized as follows. We start with a brief 
theoretical background of SQUID oscillators in Sec. ini 
In Sec. imi we describe the fabrication of a lumped el¬ 
ement SQUID oscillator and experimental measurement 
scheme. Measurements in the frequency and time do¬ 
mains are presented in Sec. IIVI A model of the electrical 
circuit is presented in Sec. El which is then used to de¬ 
rive the equations of motion of the full system. These 
equations are then reduced, accounting for the physical 
parameters used, and after incorporating thermal effects, 
are solved numerically to obtain the free evolution of the 
system for different conditions corresponding to the ex¬ 
periment. The simulated ringdowns in the time domain 
as a function of flux bias and pulse amplitude are in good 
agreement with the observations, as is shown in Sec. IVII 
Finally, in Sec. Em a scheme for using ringdown oscil¬ 
lations of a SQUID oscillator to read out a flux qubit 
is briefly discussed, followed by our conclusions in Sec. 

Eim 
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II. THEORETICAL BACKGROUND 


At the heart of all these investigations is the Joseph- 
son junction, which behaves as a nonlinear LC-oscillator 
characterized by the plasma resonance I, OJp = 
^/2^TIo/^oCJ , where, Iq is the critical current of the 
Josephson junction in parallel with its self-capacitance 
Cj and $0 = /i/2e is the magnetic flux quantum. For 
typical parameters of the fabricated junctions, is 

of the order of 100 GHz. When shunted by a large exter¬ 
nal capacitance however, the resonance frequency can be 
lowered to a few GHz for ease of performing experiments 
with an oscillator for coupling to a qubit or to fabricate 
a qubit itself. 

In this paper we consider a dc SQUID, which has two 
identical junctions in parallel, symmetrically placed on a 
superconducting loop. The dynamics of such a SQUID, 
as shown in a circuit schematic of Fig.[T] can be described 
by a two-dimensional anharmonic potential C/fc?-;-, 

given by [i,i3 

= -ih'-p^r - COS cos + 4 ~ (1) 


where Ej = /o<i>o/27r is the Josephson coupling energy 
normalizing 17, Ljo = d>o/27r/o is the Josephson induc¬ 
tance of each junction, /3 = 27rJoLg/$o = Lg/Lj^ is the 
screening parameter of the SQUID with geometric loop 
inductance Lg, if, = Ib/2,lQ is the bias current normalized 
by the critical current, and fg = d>s/<I>o is the normalized 
applied flux. = {(pA + </5b)/ 2 and p- = {pA - ^b)I‘2. 
are the two independent degrees of freedom with ipA and 
ipB corresponding to the phase differences across each 
of the two Josephson junctions (see Fig. [SKa) and Ap¬ 
pendix 1^ for a more complete description of the SQUID 
potential energy and the rest of the circuit). The sum 
of the phases across the junctions (external mode), tp+, 
couples to the current through the SQUID, while the dif¬ 
ference of the phases (internal mode), p -, couples to the 
magnetic flux applied to the SQUID. The oscillator can 
be resonantly excited by applying a short alternating cur¬ 
rent pulse to the SQUID. The pulse perturbs the poten¬ 
tial minimum of the SQUID, giving rise to oscillations of 
the phase particle about the minimum that decay at the 
characteristic frequency of the oscillator. The ringdown 
motion is mediated via the external mode of the oscil¬ 
lator and can be detected as a voltage oscillation across 
the SQUID. 


III. EXPERIMENTAL SETUP 

We investigate a lumped-element microwave oscillator 
circuit consisting of a dc SQUID shunted by a capacitor 
formed from superconducting layers. A circuit schematic 
and an optical micrograph of the device are shown in 
Fig.[TJ A microwave feedline with on-chip capacitors cou¬ 
ples signals into and out of the oscillator. Adjusting the 


bias flux <I>s as shown in the figure modulates the Joseph¬ 
son inductance of the SQUID, thus varying the resonance 
frequency of the SQUID oscillator. 



FIG. 1. (Color online) (a) Circuit schematic showing the in¬ 
put/output coupling capacitors and the SQUID oscillator, (b) 
Optical image of the fabricated circuit with a zoomed-in view 
of the SQUID with on-chip flux lines. 

Our devices are fabricated in a 5-layer process on an 
oxidized Si wafer. The initial four layers are patterned 
using photolithography, while the final layer, consisting 
of the SQUID junctions, is patterned by electron-beam 
lithography. The ground plane is formed from a 120 nm- 
thick A1 layer. The dielectric layer on top of the ground 
plane is a 150 nm-thick Si02 film deposited by plasma en¬ 
hanced chemical vapor deposition (PECVD). The Si02 
forms the dielectric for the parallel plate shunting ca¬ 
pacitor Ct and the output coupling capacitor Gout- The 
input coupling capacitor Cin (Fig.[T]) is interdigitated and 
is formed along with the microwave feedline and the top 
layer of the parallel plate capacitors in a 200 nm-thick 
A1 film. The estimated value of Gin, based on the fabri¬ 
cated finger dimensions and using the effective dielectric 
constant of a microstrip line in the standard expression 
for the interdigitated capacitor [Sll, [s^, is 0.15 pF. The 
parallel plate capacitors Gout and Ct are designed to be 
1.3 pF and 6.8 pF, respectively, based on geometry and 
film parameters, although there could be significant vari¬ 
ation because of our uncertainty in the dielectric proper- 
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ties of the Si02 film. Vias are etched through the Si02 
layer so that the bias lines and one end of the SQUID 
contact the ground plane. The SQUID loop is formed 
so that it is coupled symmetrically between the top and 
bottom plates of the shunt capacitor [Fig.[TJb)]. The ge¬ 
ometric inductanace Lg of the SQUID is calculated to be 
43 pH, from FastHenry simulations of the loop dimension 
(18 X 18^m2). 

The SQUID junctions are Al-AlO^-Al, formed by 
a standard double-angle shadow-evaporation method 
in a dedicated electron-beam evaporation cham¬ 
ber equipped with in situ Ar ion milling to ensure su¬ 
perconducting contacts between the SQUID layer and 
the junctions. The junctions are sub-micron in size, 
530 X 160 nm^, with a junction capacitance estimated to 
be 10 fF [ 3 ^. The critical current of each junction was 
estimated to be 0.4 ^A by measuring the normal state 
resistance of a nominally identical junction (684 D) and 
based on our previous characterizations of similar sized 
junctions [3^. The resonance frequency of the oscillator 
was designed to be at 3 GHz. 

Measurements are performed in a ^He cryostat with a 
base temperature of 300 mK. A schematic of the mea¬ 
surement setup is shown in Fig. [IJa). A vector net¬ 
work analyzer is used to characterize the resonance fre¬ 
quency, while the transient dynamics are studied using 
a custom-built GHz DAG (digital to analog converter) 
and a 20 GHz sampling oscilloscope. The sample chip is 
wire-bonded to a Cu stripline microwave board and en¬ 
closed in an aluminum box for magnetic shielding, which 
is anchored to the cold plate of the ^He cryostat. The 
drive-line to the SQUID oscillator is a lossy stainless 
steel semirigid coaxial cable with attenuators heat sunk 
at various stages of the cryostat to minimize noise from 
room temperature. The transmitted signals at the out¬ 
put are amplified by two High Electron Mobility Tran¬ 
sistor (HEMT) amplifiers: one at the 4K stage of the 
cryostat and another at room temperature with a com¬ 
bined gain of 70 dB. A 6dB attenuator is used at the 
output of the oscillator for 50 D impedance matching to 
the input of the cryogenic HEMT amplifier. The dc bias¬ 
ing lines have copper powder filters anchored at the 1 K 
stage of the cryostat. A cryogenic /r-metal can surrounds 
the vacuum can of the cryostat to shield the SQUID from 
external magnetic fields. 

Generation of the microwave pulses is achieved by em¬ 
ploying a EPGA based DAG board, which has been built 
based on the designs from UGSB [^. A schematic of the 
pulse generation setup is shown in Eig.j^Kb). A nanosec¬ 
ond waveform is generated digitally from the FPGA con¬ 
trolled from a computer. The waveform is passed through 
Gaussian filters and attenuators before being mixed with 
a resonant carrier tone at the IQ-mixer, producing a short 
microwave burst at the output port of the mixer. Dif¬ 
ferent amplitudes of the microwave pulse are achieved 
by varying the attenuation after the output stage of the 
mixer with a step attenuator. A typical room temper¬ 
ature trace of the microwave bursts generated by this 


setup and corresponding time domain measurement is 
discussed in the next section. 




FIG. 2. (Color online) (a) Schematic of general measurement 
setup. The network analyzer is used to measure the frequency 
response while the pulse generation setup (expanded in (b)) 
is used to measure the ringdowns in the time domain. The 
two types of measurements are separate and, as indicated by 
dashed lines in (a), the network analyzer is never connected 
simultaneously with the pulse generation setup. The SQUID 
oscillator device shown as a blue box is displayed in Fig. 1. 


IV. MEASUREMENTS OF THE SQUID 
OSCILLATOR 

The resonance frequency of the SQUID oscillator was 
characterized as a function of the applied flux at zero bias 
current (If, = 0). A network analyzer supplied a weak 
signal with a typical power of —125 dBm to the input 
of the SQUID oscillator and the 2-port complex trans¬ 
mission parameter, S 2 i{f), was measured. Although, we 
measured the complex quantity, S' 21 , the plots that we 
present only show the magnitude, |S' 2 i|- Multiple |S' 2 i| 
traces were recorded while stepping through the flux ap¬ 
plied to the SQUID. Figure [3] shows a density plot of 
the flux-modulated resonance frequency, periodic in $ 0 , 
with the color scaling (darker blue) region indicating the 
highest magnitude of |S' 2 i| and the flux axis scaled in 
units of 4>o. A phenomenological function of the form 
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fiV^s) = cos (&V$s + c) I, where V$s is the flux bias 
voltage, is used to fit the flux modulation of the resonance 
peaks over two periods. The fit parameters a, b, c are 
used to scale the flux axis and fit the frequency values, 
so that the fit curve (dashed gray line) can be plotted on 
top of the scaled resonance modulation data, as shown 
in Fig. [3l The highest frequency at integer $o occurs 
at 3.2 GHz. As expected, for /3 ^ 1, the critical current, 
and hence the resonance frequency modulate to near zero. 
The observed quality factor of the resonance is ~ 7, con¬ 
sistent with the estimated value for the output coupling 
capacitor and the chip parameters as discussed in the pre¬ 
vious Section imi In the following section, we describe a 
circuit model used to fit the measured |5'2i| traces as a 
function of applied flux by calculating the flux-dependent 
Josephson inductance based on the SQUID potential. 


A. Transmission through a SQUID oscillator with 
asymmetric coupling 

A full circuit model is shown in Fig. [SKa), where a 
Josephson junction, with parallel capacitance Cj, is in 
series with a geometric inductance Lg/2, symmetrically 
positioned on each side of the SQUID. The external shunt 
capacitance is labeled as Ct and the internal dissipation 
in the SQUID oscillator due, for example, to losses in the 
shunt capacitor, is depicted as Rt- C'in(C'out) is the input 
(output) coupling capacitor and Rz is the characteristic 
impedance at the input and output. Neglecting the ef¬ 
fects of the small junction capacitances Cj and noting 
that the dc bias current R = 0, at low-amplitude drive, 
such that the Josephson junction is nearly linear, we can 
describe the SQUID as an effective inductance 


Lt 


Ljo ^ 

2 cos 4 ■ 


( 2 ) 


Here represents the steady state value of (/?_, which 
for a given value of the applied flux fg, we can calculate 
numerically by minimizing the potential energy described 
by Eq. ([T]). 

Using Eq. m we can derive a general expression for 521 
for the effective parallel LCR tank circuit with asymmet¬ 
ric input and output coupling capacitors. A similar anal¬ 
ysis with symmetric coupling capacitors was performed 
in Ref. [ssj. We define ^21 = 2I4ut/kin such that a 
matched load of Rz = 50 0 corresponds to full trans¬ 
mission, 521 = 1 if the input were connected directly to 
the output. Then, taking the notation m || n to repre¬ 
sent the parallel impedance between impedance elements 
m and n, we arrive at 


52i(a;, Lt) 


2V),ut _ 2 II ^out) 

Vin Zout {Zt II Zont) + -^in ’ 


(3) 


where Zt is the impedance of the parallel LCR tank cir¬ 


cuit. 


^t = 




iujCi 


(4) 


and Zin = (^Rz + ^, ^out = {^Rz + 

input and output impedances respectively. 

Equation ([3]) was used to fit the measured 52 1 traces at 
each flux bias value. First, the measured 52i traces were 
scaled by the 52i transmission at low temperature that 
was measured on a separate cooldown using a direct coax¬ 
ial connection in place of the oscillator chip. This scaling 
of measured 52i of the SQUID oscillator resonance by 
the effective low-temperature baseline takes into account 
the temperature dependence in the system transmission. 

We fit the 52i(/) data by fixing /g = 0.4/rA, Cin = 
0.15pF and /3 = 0.05 and varying Ct and Rt- The esti¬ 
mates for the fixed parameters were explained earlier in 
Section iHll Measured 52i curves between ±0.3 4>o were fit 
simultaneously with the same fixed parameters and the 
best fit parameters extracted for Ct and Rt were 5.1 pF 
and 265 Q respectively. Ct is in reasonable agreement 
with our estimate based on fabrication parameters de¬ 
scribed in Section iHll We can extract a loss tangent tan 5 
for Si 02 from the equivalent resistance Rt at a frequency 
of 2.5 GHz as XjioCtRt 4 x 10“^. We note that this 
is roughly an order of magnitude larger than what was 
measured in Ref. [ 3 ^ for PECVD-deposited Si02 and 
this difference could be due to variations in deposition 
conditions and film quality. Also, a precise determina¬ 
tion of the internal loss in our capacitor is difficult with 
these measurements because the large output coupling 
strength limits the total oscillator quality factor in our 
circuit. The resonance peaks from the circuit model fits 
are plotted in Fig. [3] as solid magenta symbols. The ex¬ 
tracted ht parameters Ct and Rt are further used in the 
time-domain simulations described in Section E 

Time-domain measurements were performed by apply¬ 
ing a short microwave burst to the SQUID oscillator that 
was flux biased at 4)s = 0.3 $ 0 , where the resonance fre¬ 
quency, and hence the ringdown waveforms, have high 
flux sensitivity. An example microwave pulse generated 
by the GHz DAG is shown in Fig. UDa) along with the 
corresponding theoretical fit used for simulations of the 
ringdowns described in Section E The voltage signal at 
the output of the SQUID oscillator is shown in Fig. SDb) 
for the case when the SQUID is flux biased at 0.3 $ 0 , 
corresponding to a resonance frequency of 2.4 GHz. A 
single experimental ringdown trace shown in the figure is 
an average of 1000 traces on the sampling oscilloscope. 
The carrier frequency of the microwave burst is 2.4 GHz, 
to be on resonance when the SQUID oscillator is biased 
at 0.3 4)0 

The SQUID ringdown response to various amplitudes 
of the drive pulse was mapped out by varying the at¬ 
tenuation on the drive line of the SQUID oscillator. The 
attenuation was varied in steps of 1 dB from 0 to 40 dB at 
the top of the cryostat. At a setting of 0 dB, that is, with 
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1^211 (dB) 



(3>o) 


FIG. 3. (Color online) Density plot of |52i| vs. flux and fre¬ 
quency of the SQUID oscillator as measured from a network 
analyzer at —125 dBm power at the input of the SQUID oscil¬ 
lator chip at 300 inK. The dashed line and magenta symbols 
are from fits to the SQUID modulation as described in the 
text. The marker in black indicates the bias point where the 
pulsed measurements were taken. 


no extra attenuation other than the fixed attenuation in¬ 
side the cryostat, the measured peak amplitude of the 
burst at the top of the cryostat was 50 mV. Figure SKb) 
displays experimental and theoretical ringdown traces, to 
be described in more detail in the next section, from the 
SQUID oscillator when the pulse amplitude driving the 
SQUID is much smaller than its critical current. The 
simulated ringdown is shown as a black line while the 
experimental data is in solid circles. 


V. MODEL AND SIMULATIONS 

In Sec. El where the calibration procedure was de¬ 
scribed, and where only small drive amplitudes were con¬ 
sidered, we treated the Josephson junctions as linear el¬ 
ements, and hence the SQUID as a simple effective in¬ 
ductance. In this section we consider a more complete 
picture of the full circuit, shown in Fig. Oja), which ac¬ 
counts for the nonlinearity that is relevant when the ap¬ 
plied pulse amplitude is high. In Appendix we first 
write down the equations of motion of the full system 
in terms of five degrees of freedom. Since in our case, 
the shunt capacitance Ct is much larger than the junc¬ 
tion capacitances Cj, and the Josephson inductance of 
the junctions Ljq is in turn much larger than the ge¬ 
ometric inductance of the SQUID, Lg (a condition we 
can write as /3 = Lg/Ljg <§; 1), we can eliminate the 
fastest degrees of freedom corresponding to nodes 4 and 



t(ns) 



t (ns) 


FIG. 4. (Golor online) (a) Input pulse as measured at the 
top of the cryostat. The red dots represent the measured 
data, while the blue curve is the fit used as the input voltage 
in the simulations, (b) Example measured output voltage 
ringdown from the SQUID oscillator. Once again, the red 
dots were obtained directly from experimental data, while the 
gray curves in the backgrouird are a subset of the realizations 
obtained from the simulatioirs of the circuit. The average of 
these realizations is shown in black. 


5 in Fig. [5(a). Doing so, leaves us with a set of differen¬ 
tial equations that govern the behavior of our circuit, of 
only three degrees of freedom that treat the SQUID po¬ 
tential energy as one-dimensional. Furthermore, in Ap¬ 
pendix |Al we also discuss how, u sing the thermodynamic 
dissipation-fluctuation theorem |39j, we can include the 
effects of noise due to the non-zero temperature of the 
system. This stochastic noise turns out to play a crucial 
role in reproducing the experimental behavior when high 
amplitude input pulses are considered. Combining all of 
these factors results in the effective equations of motion 
(in units of current), which can be written using vector 
notation as 

(jigCip 4—Nn -|- /(ir = 0. (5) 

90 

Here 0 = {(pi,(p2,‘P3)’^, = [d/dipi^d/dip 2 ,d/dip^Y, 

hr = (-Un/I?z,0,0)^, n = (ni, 712,71-3)^ and (j)o = 
$o/27r. Each variable ipi^ represents the superconducting 
phase at node I, with the corresponding voltage defined 
as (j)Q0i. The thermal noise in the circuit is modeled by 

including a current noise source of strength J in 
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(a) 




y^Ui = P ( (sin(47r/s + 2(^2) + siii(2(/j2)) | • (H) 


Eq. ([5]) forms a set of stochastic differential equations 
that we can numerically solve for any although each 
solution only gives us a single realization. Averaging over 
many such realizations (500 in our case) produces a curve 
that can be directly compared to the experimental data, 
which we do in Sec. m In our simulations we chose the 
initial conditions that correspond to the system being at 
rest, near the potential energy minimum, which we can 
express as (^(0) = 0 and (p(0) = (0, —/sTt, 0)^ respec¬ 
tively. We then let the system thermalize by evolving 
Eq. [5] without an external drive present (/dr = 0), with 
only the thermal noise influencing the evolution. In the 
final step, the input pulse shown in Fig. UKa) is applied, 
which in turn excites the system, leading to ringdown 
oscillations. 


A. Model Limitations 


FIG. 5. (Color online) (a) Full circuit model of the SQUID 
oscillator of Fig. [T]and (b) reduced circuit, valid in the limit 
of small amplitudes when the Josephson junction responds 
linearly. In this regime the SQUID is treated as an effective, 
flux-dependent inductance Lt- In both (a) and (b), the blue 
numbers represent node labels. 

parallel with each resistor i?^. We take ks to represent 
the Boltzmann constant, T the temperature and each n^, 
a normally distributed random variable, which satishes 

(n^(t)) =0 (6) 

(nPt)ny(t')) (7) 

In our case, taking Cs = Ct + Qn + Cout, the matrices 
corresponding to C and can be written as 

/ Qn -Qn 0 \ 0 0\ 

c = -Qn Cs -Cout ,R-^ = 0 ^ 0 , 

V 0 -Cout Cout / V 0 0 

^ ( 8 ) 

and N as simply 

N = ^/2kBTR-l. (9) 

Finally, UeS represents the effective (drive free) poten¬ 
tial energy of our system, which can be decomposed as 
Ues = Uo + Ui- The terms Uq and Ui represent the con¬ 
tributions, up to zeroth and first order in /3 respectively 
(as explained in Appendix , which result in 

V^Co = I 2£’jcos(7r/s)sm(v32+ 7r/s) 1 , (10) 


In our model, we neglect the internal dissipation as¬ 
sociated with each Josephson junction. This is reason¬ 
able when the amplitude of the excitation applied to the 
junctions is smaller than their critical current, as the ef¬ 
fective resistance shunting each junction in this case is 
large enough that its effect on the damping of the junc¬ 
tion phase can be neglected. However, in the instances 
when the driving current exceeds the critical current, 
each junction experiences a resistance that can be of the 
order of its normal-state resistance Rn [1^ , which for the 
SQUID oscillator studied here is 68417. Nevertheless, in 
our circuit, the dominant source of noise is the 5017 out¬ 
side load that couples to the SQUID oscillator via Cout, 
as can be seen from Fig. [SJ When this load is mathe¬ 
matically transformed as impedance parallel with the os¬ 
cillator over the frequency ranges of the input pulses we 
apply, its resistive component is never more than 10017, 
hence a few times smaller than all other sources of noise 
in the system such that its effect is by far the most dom¬ 
inant. 

Furthermore, we neglect any quantum corrections to 
the noise correlation function and, as is shown in Eq. o, 
treat it, just as the rest of the system, fully classically. 
This is typically a reasonable assumption in the limit of 
fujj <C 2kBT with ui being the applied, flux dependent, 
effective natural frequency of the oscillator circuit. In the 
case of the experimental parameters used here, this limit 
is largely satisfied, although in the worst case, when the 
flux through the SQUID is close to integer multiples of 
a flux quantum (where the effective natural frequency of 
the oscillator is largest) we are slowly approaching a case 
where hu < 2kBT^ in particular with T = 0.300 K and 
at /s = 0, we have 2 ^t i- Calculating the leading 
order correction to the quantum version of the correlation 
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function [iO, EH , leads to the following adjusted Eq. dZD 

{n,{t)nj{t')) = 6{t - ( 12 ) 

which in turn we can interpret as a factor of ^ 1/24 
change in the amplitude of the noisy current sources as¬ 
sociated with each of the resistors. Given that this is the 
worst case, and this correction gets smaller as the flux 
that is threaded through the SQUID shifts away from 
integer multiple of $o, we neglect it in our simulations. 


VI. VOLTAGE RINGDOWNS 

While we can simulate the evolution of an arbitrary 
degree of freedom, of particular interest is $3 = as 

it corresponds to the output voltage, which is precisely 
what is measured in the experiments. As already dis¬ 
cussed, we have two “knobs” that can be controlled in a 
given experimental run; the flux bias /s, and the ampli¬ 
tude of the input pulse. The rest of the parameters are 
fixed at the values described in Sec. El To explore the 
behavior of our system better, it is therefore instructive 
to vary one of these control knobs while keeping the other 
constant. 


A. Amplitude Scans 


V (mV) 




4 

2 

0 



V(mV) 

1 ^ 

■■ 2 

■■ 0 



We first look at amplitude scans, where we fix the 
applied flux bias fs and vary the amplitude of the mi¬ 
crowave burst, recording a ringdown trace for each excita¬ 
tion strength of the SQUID oscillator. Figure H] displays 
density plots of such a case, obtained with experimen¬ 
tal data (top plot) and from simulations (bottom plot). 
The flux fs is fixed at 0.30, which corresponds to the 
natural frequency of the SQUID of 2.4 GHz (when the 
SQUID is operated in a linear regime), satisfying a res¬ 
onance condition. The y-axis, is the attenuation setting 
from the highest (40 dB) to the lowest attenuation (0 dB) 
corresponding to increasing burst amplitude towards the 
bottom of the plot and the x-axis is the ringdown time 
in nanoseconds. The color scale indicates the amplitude 
of the ringdown. With increasing pulse strength, the fre¬ 
quency of the ringdowns decreases. This becomes partic¬ 
ularly pronounced for attenuation levels less than 20 dB. 
This shift to lower frequencies arises because the fictitious 
particle, whose position coordinate can be described by 
the (f )2 degree of freedom, begins to explore the nonlinear 
(flatter) part of the potential energy landscape. For at¬ 
tenuation levels less than 12 dB, we observe a sharp drop 
in the resulting ringdowns. At this point, the strength 
of the drive is now of the order of the critical current of 
the SQUID. The stochastic nature of the thermal noise 
causes different realizations to escape the potential well 
at different times, which in turn causes a substantial de¬ 
cay in the ringdown signal strength. This phenomenon 
is discussed in more detail in the next section. 


FIG. 6. (Color online) A comparison of the amplitude depen¬ 
dence of ringdowns at a non-integer flux bias of /s = 0.30, 
for the input signal of Fig. 14(a) at 2.4 GHz. The plot in (a) 
shows experimental data, while in (b), the corresponding sim¬ 
ulations. The attenuation on the drive pulse in dB is shown 
on the y-axis, with decreasing numbers implying increasing 
amplitude of the input pulse. The ringdown time is shown 
on the x-axis while the ringdown amplitude is represented 
by the color scale. The purple dashed line in (a) indicates 
the amplitude corresponding to the ringdown trace shown in 
Fig. mb). 


B. Escape from the Potential Well 

The problem of particle escape from a potential well 
due to thermal noise has been investigated both theo- 
reticall y a s well as experimentally in a variety of stud¬ 
ies [42l - l4^ . In the case of a dc SQUID, this rate can 
be approximated to be proportional to D exp (—Ub/ksT) 
where Ub represents the potential energy barrier height 
that the particle has to overcome, and D the natural fre¬ 
quency along the direction of escape. In our case, since 
we do not “tilt” the potential with a dc bias current, 
the escape time (inverse rate) can be shown to be much 
larger than the typical experimental run time. This is 
true over almost all settings of the applied flux bias fs, 
except when fs ~ 0.50, where the potential barrier is 
close to being flat. 
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FIG. 7. (Color online) Low and high drive response. Each row shows plots of 4>o<p2{t) and (f>oipi(t) (the output voltage) 

respectively. In plots (a-c), the amplitude of the input pulse is low with respect to the critical current of the Josephson junctions, 
whereas in (d-f) it is high. The gray curves show a small subset of individual realizations used to calculate the averages (blue 
curves). In the case of the third column we also show the experimental data for the same parameters (red dots) and how it 
compares to the simulation results. From the plots, one clearly sees how when the drive is low, all the realizations stay within 
the same potential well (see plot a), and their relative phase shift only varies slightly. In the case of high drive, we observe that 
different realizations tend to end up in different wells at different times (as shown in plot d), which introduces a relative phase 
shift between them. This in turn leads to faster decay of the average voltage, which is reflected in plots (e) and (f). 


Thermal fluctuations, however, still end up playing an 
important role in the evolution of the system. In par¬ 
ticular, we find that during strong pulses that excite the 
system to amplitudes in the vicinity of the SQUID’s crit¬ 
ical current, the thermal fluctuations can cause a strong 
mixing in the phases of various realizations, resulting in 
a damping of the ringdowns. To our knowledge, no de¬ 
tailed analytical study of this effect, with strongly time- 
dependent, transient pulses has been performed yet, but 
we can still study the situation numerically. To do this, 
we once again fix the applied flux bias at /s = 0.30 as 
in the previous section, and concentrate on two differ¬ 
ent pulses: the first at an attenuation of 15 dB and the 
other at the attenuation of 10 dB. From Fig. [51 we can 
see that these correspond to cases where substantial ring- 
down voltage is observed (the former case) and where the 
ringdowns are dramatically suppressed (the latter case). 

In order to understand this behavior in more detail, we 
look at the evolution of the individual realizations, that 
so far have been averaged to obtain results comparable 
with the experiment. We stress, however, that while the 
behavior of the full circuit is largely governed by the dy¬ 
namics of the SQUID, the experiment only provides us 
access to the external voltage — the voltage at node 3 in 


Fig. O — which in our simulations is represented mathe¬ 
matically as 4>o<f3. To directly observe the stochastic na¬ 
ture of the escape from the potential well, we need to look 
at the individual realizations of the full system. Of par¬ 
ticular interest are the following simulation variables: ip 2 , 
which represents the phase (i.e., the effective “position”) 
of the SQUID degree of freedom, that dominates the evo¬ 
lution of the system, (f>o'f 2 , which represents the voltage 
across the SQUID (or alternatively an effective “velocity” 
of the particle in the well), and finally which is the 

voltage that we can directly compare to the experimental 
data. Figure [7] shows plots that describe the evolution of 
these variables as a function of time. The top row shows 
data for a case of the low amplitude, 15 dB attenuation 
pulse, while the bottom row shows the case of high ampli¬ 
tude, 10 dB attenuation pulse. The leftmost column rep¬ 
resents the middle column (^o<P 2 (t), and finally the 

rightmost column is the output voltage, namely (poip 3 {t). 
In each case, the grey curves represent a subset of re¬ 
alizations that are averaged (curves in blue). The red 
dots in the plots from the rightmost column represent 
experimental data for the same set of parameters as the 
simulations. The key signature of the escape can be seen 
in the leftmost column. Here, when the pulse amplitude 
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is low (top row), virtually all the realizations stay within 
the same potential well — as one can see by noting that 
they all oscillate around the same value of ^p 2 ^ namely 
ip 2 ~ —O.SOtt. In the case of the high amplitude pulse 
(bottom row), different realizations jump out to differ¬ 
ent potential wells. The stochastic nature of the noise 
causes these jumps to happen at different times, which 
leads to a randomly shifted phase, as well as a different 
steady state value of Lp 2 ■ This has a substantial effect on 
the “velocity” (or 02 ) of these realizations, as shown in 
the central column of Fig. [T] The result is a dramatic 
randomization in the phase of 02 , and as a result, of 03 , 
which is proportional to the output voltage of the circuit. 
As we see from the experimental voltage (red dots), the 
agreement of the measured data with the simulations is 
good. Finally, we stress that including the stochastic ef¬ 
fects of thermal noise in our simulations has been crucial 
in reproducing this behavior. 


C. Flux Scans 

We explore the voltage ringdown behavior further by 
studying their dependence on the magnetic flux applied 
to the SQUID. Here, the amplitude and frequency of the 
microwave burst is fixed while we vary the flux applied to 
the SQUID through one period of a flux quantum. The 
pulse frequency once again is chosen to correspond to res¬ 
onance at the SQUID flux bias fs 0.30. The density 
plots of the flux-modulated ringdown traces is shown in 
Fig. [5] for three different pulse amplitudes. The top row 
shows plots obtained from experimental data, while the 
bottom row shows the simulations. The leftmost column 
has a low input pulse amplitude with 20 dB of attenua¬ 
tion, well below the critical current of the SQUID, the 
central column shows data for an input pulse with 15 dB 
of attenuation, while the rightmost column has a high 
amplitude pulse with 10 dB of attenuation. By varying 
the applied flux through the SQUID, we are changing 
its effective inductance, and hence its natural frequency. 
It is worth stressing that this nonlinear dependence of 
the natural frequency on the applied flux is true even 
in the limit of small oscillations of the SQUID (where 
\^P 2 \ <SC 1), as was already discussed in Sec. |IVl Let us 
first concentrate on the leftmost column of Fig. |S1 Here 
the drive amplitude is still small and the nonlinearity of 
the potential energy in ^p 2 is only beginning to play a 
role. Yet, as the applied flux bias fs varies between 0 
and 1, the ringdowns tend to fan out. As expected, the 
amplitude is largest near the flux bias of fs ~ 0.30, since 
this is where the SQUID is resonant with the input pulse, 
and it is suppressed elsewhere. The results are also con¬ 
sistent with the fact that the natural frequency, up to 
zeroth order in /?, is proportional to a/ cos(7r/s). Hence, 
near /s ^ 0, the variations in the ringdown structure are 
small, while at the same time, one sees a very abrupt 
suppression near fs ^ 0.50. Here the effective natural 
frequency of the SQUID is very small and the short in¬ 


put pulse is unable to induce strong oscillations. This 
is an adiabatic regime, where the excitation of the cir¬ 
cuit strongly follows the input pulse. In this regime, the 
ringdown suppression due to a highly off-resonant pulse 
can be confirmed further by studying individual realiza¬ 
tions and showing that they stay in the same potential 
energy well as they started in, in contrast to what is ob¬ 
served during an escape — see Sec. IVI Bl Furthermore, 
a very similar ringdown structure can be obtained in a 
case where a simple harmonic oscillator, with the same 
flux-dependent form of natural frequency, is driven with 
the same pulse waveform. The situation is largely similar 
in the middle column of Fig. [8] The key difference here 
is that now, not only is the natural frequency of the sys¬ 
tem nonlinear in the applied flux, but the amplitude of 
the input pulse is large enough for the SQUID degree of 
freedom ip 2 to start exploring the nonlinear regions of the 
potential energy well. This effect is particularly strong 
around the applied flux bias for which the SQUID is res¬ 
onant with the input pulse (near fs ^ 0.30). This in turn 
affects the degree of variation of the ringdown frequency 
with respect to fs, as can be seen in the plots. Finally, in 
the rightmost column we see a case of a strongly driven 
system. The resulting plots show an overall suppression 
of ringdown oscillations across all values of fs , when com¬ 
pared to the instances with smaller drive amplitudes. In 
this case, the reason is two-fold. By once again studying 
the individual realizations as in Sec. IVIBl we can con¬ 
clude that for the applied flux away from fs = 0.50, the 
main culprit in the suppression is the randomization of 
the phase of (p 2 due to the stochastic escape from the 
potential well. Near fs = 0.50 however, as in the case of 
low amplitude pulses, the main reason for the suppres¬ 
sion is the off-resonance condition, where the frequency 
of the pulse is much greater that the natural frequency 
of the SQUID. 


VII. APPLICATION TO FLUX 
MEASUREMENTS 

The high sensitivity of SQUIDs to applied flux has 
made them exquisite detectors of flux signals. They have 
been useful in various metrology experiments [^, 
and more recently have played an important role in the 
field of quantum computing, as measurement devices for 
flux qubits [lOl - fTsl . [^ . In these applications, a flux qubit 
is typically coupled inductively to a SQUID and hence 
affects the net applied flux that is threaded through the 
device. One of the original SQUID-based readout ap¬ 
proaches involves biasing the SQUID with an appropri¬ 
ately selected dc current, such that the SQUID is put 
in a running state with non-zero voltage, with a high 
probability if the qubit is in one state, and with negligi¬ 
ble probability when the qubit is in the other state (7|. 
In an alternate approach, the SQUID is driven with a 
continuous sinusoidal signal while monitoring the result¬ 
ing phase shift, which is qubit-state dependent [s^. Yet 
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FIG. 8. (Color online) Fixed amplitude flux scans with the applied flux bias fa between 0 and 1, and with the input signal of 
Fig. Si: a) at 2.4GHz. The top row (a-c) shows experimental data, while the bottom row (d-f), uses data obtained by running 
stochastic simulations. The amplitude of the input pulse increases from left to right, with the leftmost column showing results 
for 20 dB attenuation pulses, the middle column for 15 dB attenuation pulses, and finally the rightmost column for 10 dB 
attenuation pulses. As the amplitude increases, one clearly observes the effects of the nonlinearity of the system. See main text 
for a more detailed discussion. 


another proposed readout scheme uses brief, but strong, 
current bias pulses to the SQUID, which result in ring- 
down dynamics with an amplitude, and possibly phase, 
that depend on the qubit state [53l| . This is somewhat 
analogous to what is presented in the experiment de¬ 
scribed here, although clearly in our case the flux dif¬ 
ferences are due to a global flux biasing as no qubit is 
actually present. Furthermore, the theoretical proposal 
outlined in shows a full quantum treatment of the 
qubit and SQUID system, but only considers an ultra- 
short dc pulse, much shorter than the inverse character¬ 
istic qubit frequency. Due to the relatively high tem¬ 
perature of our measurements, our present experiment 
is in the classical regime, and involves input pulses with 
a time scale comparable to the dynamics of the circuit. 
Nevertheless, it is still useful to explore briefly just how 
the process of discrimination between two or more differ¬ 
ent flux states could be accomplished with our system. 
First, one can expect that the total applied flux through 
the SQUID would consist of some static bias flux $bias 
plus a signal flux that is to be measured, say ^signal- One 
could then send a microwave pulse through the SQUID, 


analogous to what was considered here, and record the 
corresponding ringdown voltage. Some post-processing 
of this ringdown voltage, such as, for example, taking its 
root-mean-squared value, Ums, integrated over a suitably 
chosen time range, would provide a level corresponding 
to the signal flux. As long as these levels of various val¬ 
ues of <1’signal Can be distinguished, one has an effective 
flux meter. For a given input pulse, assuming that Ums 
is a well behaved function of the total SQUID flux $s, 
one possible way to find the best ‘I’bias over a range of 
flux where Ums is monotonic would be simply to look 
for the largest slope of I4ms with respect to $8, namely 
maximizing 9Urms(iint,‘f’s)/9‘fs over all possible flux $8 
between 0 and O.S^q (due to symmetry) and integration 
times tint • This would give the greatest contrast between 
the cases of $bias + ^signal and d>bias - ^’signal- Figured] 
shows an explicit example of this kind of flux discrimi¬ 
nation based on our measurements, where we calculate 
Ums over a time range between 2.1 and 3.4 ns. The data 
that is being used corresponds to the input signal with 
20 dB attenuation, and is the same as in the flux-scans 
from Fig. jS] The dots represent results obtained from 
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FIG. 9. (Color online) Root mean square output voltage Kms 
as a function of the applied flux bias cl’s, calculated over a time 
range between 2.1 and 3.4 ns. The plot uses data obtained 
with the 20 dB attenuation input pulses, and is the same as 
in the left column of the flux-scans from Fig. [H] The dots 
represent results calculated from the experimental data (top 
row, leftmost column), while the solid line is produced using 
the simulations (bottom row, leftmost column). By biasing 
the flux through the SQUID near a point where the slope is 
high, for example at ~ 0.36'Fo, one can have a means of 
distinguishing between different flux signals — see main text 
for more details. 


the experimental data (top row, leftmost column), while 
the solid line is produced from the evolution calculated 
through stochastic simulations (bottom row, leftmost col¬ 
umn). We can further make a crude calculation of the 
required sensitivity that one would need with the data 
from Fig.|9]to distinguish between two hypothetical flux 
qubit states. Setting the bias flux at $bias ~ 0.36$o the 
slope is roughly 5mV/$o- If we assume a conservative 
noise temperature of 150 mK for a ^ 3 GHz amplifier with 
a bandwidth of 100 MHz [s^, the rms voltage noise 
at the amplifier input would be ^ 200 nV. If we take the 
5 mV/$o slope for the signal at the output extracted from 
Fig.l and divide by the net gain of the HEMT amplifiers 
(~ 55 dB), this becomes 9/iV/$o at the SQUID oscilla¬ 
tor output. We consider a peak-to-peak qubit flux signal 
of 22md>o, which is reasonable!^ [^. . considering the 

back-action on the qubit would also likely be significantly 
less compared to a switching dc SQUID measurement, 
since the SQUID never enters the running state. This 
then corresponds to a SNR of ^ I. So, we would be right 
at the threshold for reading out the ringdowns and dis¬ 
tinguishing between the two qubit states in a single shot. 
We should further stress that one could likely do better 
by both using more sensitive amplifiers and optimizing 
various parameters. Of particular importance would be 
integration time tint, the pulse amplitude, as well as the 
quality factor of the SQUID oscillator, all of which the 
Urms curves are highly dependent on. 


VIII. CONCLUSIONS 

In conclusion, we have studied the transient behavior 
of a dc SQUID operated as a nonlinear oscillator under 
pulsed ac excitation. Both experimentally as well as nu¬ 
merically, we applied signals of various amplitudes for 
different flux bias, while observing the resulting voltage 
ringdowns. In order to account for the non-zero temper¬ 
ature of the experiment, we used the Johnson-Nyquist 
approach and modeled resistors as noisy current sources. 
This let us numerically reproduce the stochastic escape 
dynamics observed when the SQUID was driven with 
high-amplitude pulses. Finally, we briefly discussed the 
potential applicability of our system, and in particular 
the observed ringdown dynamics, to flux measurement. 
We found a good general agreement between the experi¬ 
mental data and results obtained through numerical sim¬ 
ulations. 
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Appendix A: Equations Of Motion 

In this section we present a derivation of the equations 
of motion of a circuit which was used to model our 
experimental apparatus. We start with a full description 
of the system at zero-temperature and reduce the equa¬ 
tions of motion by eliminating fast degrees of freedom. 
We then add the effects of the temperature-dependent 
noise. 


1. Zero Temperature 

As discussed in Section imi above, the circuit diagram 
is shown in Fig. [SJ Our model assumes that the exter¬ 
nal flux is delivered directly to the SQUID loop, and 
other branches have no intrinsic geometric inductance. 
We further neglect the mutual inductance in the system 
other than the one that mediates the external flux $s. 
To obtain the equations of motion, we follow the treat¬ 
ment of Devoret [s^. With each node f, we associate a 
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corresponding node flux <i?i related to a node voltage by 
We express the currents across ele¬ 
ments in terms of and using Kirchoff’s current conser¬ 
vation conditions at each node i, arrive at the equations 
of motion 

=C'in($i -$ 2 ) 

.... 2 2 

- $ 2 ) =—($2 - d>4 + $s) + ^(^>2 - 
-^9 ^9 

+ Ct^2 + ■5“‘i’2 + C'out(^2 ~ ^ 3 ) 
C'out(^2 - ^ 3 ) =^^3 

tiz 

2 1 . .. 

— ($2 - «>4 + =h sin($427r/$o) + ^^>4 + C'jd>4 

Ijg iCi 

2 1 • •• 

($2 — $ 5 ) =/q sin($527r/$o) -|- + C'j^s 

-L/g lli 


Next, taking the flux quantum <i>o = 2.07 x 10“^^Wb = 
we perform a change of variables so that = 
Here, a difference (pi — ipj for some i j, 
corresponds to the superconducting phase difference. We 
further take <p± = 5 ( 7’4 ± T’s), C'e = C'in + (llout + Ct, 
Ljo = d>o/27r/o and rewrite the external flux in terms 
of the ratio /s = After dividing all equations by (po, 
we have 


0 =Cin‘Pl - Cin<fi2 + \ Vin 

Hz 00 Az 

0 = — Cin^il -|- C'£,(p2 — Cout03 + 

rit 


+ y- (7’2 - </?+ + tt/s) 

0 = — Cout'^2 + C'out'^3 + 7^03 

tiz 

2 4 

0 =2Cj(pj^ -V -— sin (^+ cos (p_ - — - <f++ T^fs) 


Ljo 


L„ 


2 4 

0 =2Cj(p- + -— sin(p_ cos(/9+ - — {-p- + irfs). 
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Thus, we end up with equations of motion for five de¬ 
grees of freedom. In order to further simplify the above, 
we note that in our case, the capacitances (or effec¬ 
tive masses) of oscillators (/?+ and (p- are two orders 
of magnitude smaller than that of ip 2 - Furthermore, 
the Josephson inductance Ljq is much greater than the 
geometric inductance Lg. This allows us to apply an 
Oppenheimer-Born-like approximation and eliminate the 
fast-oscillating degrees of freedom ip+ and tp-. To do 
this, we first define a potential energy U that can be as¬ 
sociated with our system. Neglecting terms due to the 
external drive, we have 


U 


- cos (p+ cos ^- + -p if - - T^fsf 

+ ^ - (7’2 + tt/s))^ , 


(Al) 


with P = Lg!L jo- Next, we fix the slow variable (p 2 , and 
note that since /3 <IC 1, the second and third terms in 
Eq. dSU will dominate. Hence the minima of U will be 
close to (p+ = tp 2 + T^fs and tp- = tt/^. By expanding U 
near these points and minimizing, we can calculate the 
corrections to the minimum points. Keeping terms up to 
first order in /3 we arrive at 

*>2'" = nf. - (A2) 

= + (A3) 


These results are then substituted back into the expanded 
potential energy, which leads to UeS = Uo + Ui, now only 
in terms of p 2 and with 


Up 

2Ej 


cos (tt/a) cos {P 2 + tt/s) 


and 


^ = - ^ (sin^ (tt/s) cos^(7r/A -F P2) 

+ cos^ (tt/a) sin^(7r/s + P 2 )) ■ (A4) 


We have distinguished between the contributions to the 
effective potential energy between terms of different or¬ 
ders in /3. Uo neglects the geometric inductance com¬ 
pletely, while Ui shows the correction up to first order 
in /3. We can hence write a set of effective equations 
of motion with UeS as the potential energy, while now 
including drive term, as 


0 =CinPl - CinP2 + 1^01 - , Mn 

Ez (poEz 

0 = — CinPl + Cs‘P2 — CoutPS + 75“7’2 

Et 


-‘JO 


-13 


■cos (7r/A)sin(((92 -|-7r/A) 
sin(47r/A -|- 2p2) + sin(2(/92) 


2L 


JO 


0 — — CoutP2 + CoutPO + 7^03 ■ 

tiz 


(AS) 


(A6) 

(A7) 


The next step is to account for the non-zero temperature 
of the system. 


2. Non-zero Temperature 

We find that the accounting for thermal noise is of par¬ 
ticular importance when comparing with the behavior of 
the experimental system in our simulations, in particular 
at high amplitude pulses. In order to model these ef¬ 
fects, we use the thermodynamic dissipation-fluctuation 
relation [s^. Thermal noise in the circuit is modeled by 

including a current noise source of strength J in 
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parallel with each resistor Ri. We take ks as the Boltz- 
man constant, and T the temperature of the system. Fur¬ 
thermore, each ni{t) represents a normally distributed 
random variable, namely rii € A/’(0,1), that satisfies the 
following 


Adding such a noisy current in parallel with each of the 
resistors to the effective model derived in Sec. ixn leads 
to classical Langevin equations, which after rearranging, 
can be written in a vector form as Eq. ®. 


{ni{t)) =0 (A8) 

{ni{t)nj{t')) =S{t -t')5ij. (A9) 
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